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In this paper, we consider a continuous-time autoregressive fractionally integrated moving av- 
erage (CARFIMA) model, which is defined as the stationary solution of a stochastic differential 
equation driven by a standard fractional Brownian motion. Like the discrete-time ARFIMA 
model, the CARFIMA model is useful for studying time series with short memory, long mem- 
ory and antipersistence. We investigate the stationarity of the model and derive its covariance 
structure. In addition, we derive the spectral density function of a stationary CARFIMA process. 
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1. Introduction 

It is well known that the long-range dependence properties of time scries data have diverse 
applications in many fields, including hydrology, economics and telecommunications (see 
Bloomfield [8], Sowell [35], Beran [4], Robinson [32], Baillie [3] and Ray and Tsay [31]). 
Autoregressive fractionally integrated moving average (ARFIMA) models are a popular 
class of discrete-time long memory processes (see Granger and Joyeux [19] and Hosking 
[21]). For continuous-time long memory modeling, see Viano et al. [40], Chambers [13], 
Comte [14], Comte and Renault [15], Brockwell and Marquardt [12], Tsai and Chan [37] 
and Tsai and Chan [38]. 

In contrast to long-range dependence, which corresponds to the singularity of the 
spectrum at the origin, antipersistent time series are covariance-stationary processes with 
zero spectral density at the origin (Dittmann and Granger [16], Bondon and Palma [9]). 
Interesting examples of antipersistent processes include: Kolmogorov [24], which models 
the local structure of turbulence in incompressible viscous fluids; Ausloos and Ivanova [2], 
which models the temporal correlations of fluctuations in the Southern Oscillation index 
(SOI) signal; Simonsen [34], which measures correlations in the Nordic electricity spot 
market. For financial applications of antipersistent processes, sec Peters [29] and Shiryaev 
[33] , both of which model the implied and realized volatility of the S & P 500 index. For 
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further examples of antipersistence modeling, see Otway [28], Beran and Mazzola [7] and 
Metzler [26]. 

Beran and Feng [5] observed that antipersistent processes can be generated by overdif- 
ferencing. In other words, a slightly non-stationary process lying between a stationary 
long-memory process and a random walk becomes antipersistent after the first differenc- 
ing. For example, Beran et al. [6] found that the first-differenced daily world copper price 
from January 2, 1997 to September 2, 1998 comprised an antipersistent process and a 
non-zero mean function. This implies that the original process was a trend function, plus 
a nonstationary error process whose first difference was antipersistent. Karuppiah and 
Los [23] also noted that many intra-day foreign exchange rate series are antipersistent. 

Recently, Tsai and Chan [37] proposed the continuous-time autoregressivc fractionally 
integrated moving average (CARFIMA) model, which is defined as the stationary solu- 
tion of a stochastic differential equation driven by standard fractional Brownian motion 
with the Hurst parameter 1/2 < _ff < 1. The model is useful for studying time series data 
that exhibits long-range dependence properties. In [37], Tsai and Chan derived the au- 
tocovariancc function of the stationary CARFIMA process and in [38], they derived its 
spectral density function; both functions are derived under the condition 1/2 < _ff < 1. 

In this paper, we focus on the stationarity, autocovariancc function, spectral density 
function and other properties of the CARFIMA process with the Hurst parameter < 
H < 1. The CARFIMA process with 0<-ff<l/2isan antipersistent process that is 
useful for studying time series data exhibiting antipersistent properties. For H = 1/2, the 
CARFIMA process becomes a continuous-time autoregressivc moving average (CARMA) 
process. We present the main results in Section 2, followed by their proofs in Section 3. 
Finally, in Section 4, we make some concluding remarks. 

2. Main results 

Heuristically, a CARFIMA(p, H, q) process {Yt} is defined as the solution of a pth order 
stochastic differential equation with suitable initial conditions. It is driven by a standard 
fractional Brownian motion with the Hurst parameter H and its derivatives up to and 
including the order 0<q<p. Specifically, for t>0, 

yI"^ - a,Y}P-'^ a,Yt-ao (1) 

where {Bt^H = ,t > 0} is a standard fractional Brownian motion with the Hurst 
parameter < H < 1, dY^'' = Y^^^ dt,j = 1, . . . ,p — 1. The superscript '^^^ denotes a 
j-fold differentiation with respect to t. We assume that a > 0, ai 7^ and (3q ^ 0. 

Let < H < 1 be a fixed number. It is well known (see, e.g., Duncan et al. [17]) that 
there exists a Gaussian stochastic process {B^ ,t > 0}, which satisfies the following three 
properties: 

(a) its initial condition is Bq = 0; 



180 



H. Tsai 



(b) it is of zero mean, that is, E(i?^) = for all t>Q; 

(c) its covariance kernel is defined as 

E(i?fBf) = i{|tr+i.r-ii-.r} 

for all .s,i > 0. The Gaussian process {B^} is called a standard fractional Brownian 
motion with the Hurst parameter H . A standard fractional Brownian motion with H = 
1/2 is the same as a standard Brownian motion. In addition, a fractional Brownian 
motion defined over non-negative t can be extended and defined for all real numbers. 
Specifically, for H G (0, 1), we have, for — cx) <t < oo, 

/oc 
{{t-u)"^-'^' -{-u)'^-'^'}dBu, 
-oo 

where Dh = [2i7r(3/2 - H) /{T{H+ 1 /2)r(2 - 2H)}Y/^, T{-) is the Gamma function and 
{Bu} is a two-sided standard Brownian motion with covariance F,{BsBt) = min(|s|, \t\) if 
s and t have the same sign and if s and t have different signs (see Taqqu [36]). 

Let Z he & set of integers. The stationary process {Ft}t£Z^ defined by Ft ~ B^^ — B^ , 
is called fractional Gaussian noise. The autocovariance of the noise is 

^F{k)^^{\k + 1\'" -2\k\^" + \k-l\'"}. 

It can be easily demonstrated (see, e.g., Taqqu [36]) that for ky^O, 7F(fc) = if iJ = 1/2, 
7F(fc) < ff < < 1/2 and 7f(A:) > ff 1/2 < < 1. For H ^ 1/2, 

-/Fik)^H{2H-l)\k\^"-^ asfc^oo. 

For 1/2 < H <1, J2'k'=-oo'^F(^) ~ noise is said to be persistent, to have long 

memory or to be long-range dependent. For < < 1/2, J2T=-oo Ipi^) = 0; therefore, 
the process is negatively autocorrelated at all positive lags and the noise is said to be 
antipersistent. For H = 1/2, all correlations of the process at non-zero lags are zero, that 
is. Ft and F^ are uncorrelated for s. The spectral density of the noise is given by (see 
Beran [4]) 

oo 

/f(w) = {2(l-cosw)} ^ |w + 2/c7t|-2ff-i. 

k— — oo 

Note that fpiuj) is 0{\uj\^~^") for uj^Q. 

Taqqu [36] considers continuous- time fractional Gaussian noise {Ct^t^R, where R is 
a set of real numbers, namely Ct = -Bj+i — B^ , and shows that the spectral density of 
continuous-time fractional Gaussian noise is given by 

fc{oj) = ^a2r(2i/+ l)sin(7ti/)|e'- - llVl"'''"' 
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Fractional Brownian motion is not difFerentiable (Mandelbrot and Van Ness [25]), so 
the stochastic equation (1) has to be appropriately interpreted as an integral equation, 
as explained below. Analogous to continuous-time ARMA processes (sec, e.g., Brockwell 
[10]), equation (1) can be equivalently cast in terms of the observation and state equations 

Yt^P'Xt, i>0, (2) 
dXt = {AXt +ao6p)dt + aSp dSf , (3) 

where the prime superscript denotes taking the transpose. 
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and Pj = for j > q. See Tsai and Chan [39] for a proof of the equivalence between 
equation (1) and the set of equations (2) and (3). The process {Yt,t > 0} is said to be 
a CARFIMA(p, H, q) process with the parameter {9, a) = {ao, . . . ,ap, (3i, . . . , f3q, H,a) 
if Yt = (3'Xt, where Xt is the solution of (3) for the initial condition Xq. Similarly to 
equation (5) of Tsai and Chan [37], the solution {Xt,t > 0} of (3) can be written as 

Xt = e^'Xo + ao f e^^'-^^Sp du + a f e^^'-^^p dSf , (4) 
Jo Jo 

where 

00 

e^*=/, + 5]{(At)"(n!)-i} 

71 = 1 

and Ip is the p x p identity matrix. The stochastic integration in (4) is defined in terms 
of the limit of Ricmann sums because it only involves deterministic integrands. 

For a random initial condition Xq, the mean vector of {A"t}, denoted by fix,t, satisfies 
the following equation: 

Aix.t = e^Vx,o + — (c^*-/p)<5i, (5) 
ai 

where 61 = [1, 0, . . . , 0]'. If fix,o is chosen to be —{ao/ai)Si, then iix,t becomes 
— {ao/ai)Si, which is independent of t. If all the eigenvalues of A have negative real 
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parts, then it can be easily shown that (2) and (4) imply that, when t — > cxd, Yt converges 
in distribution to a normal random variable with mean — ao/o^i and variance Vy, where 
Vy =7y(0) is defined in equation (11). Thus, the stationary solution, if it exists, must 
be Gaussian. 

The stationary CARFIMA process defined over non-negative t can be extended so that 
it is a stationary process over all real t. For simplicity, wc assume that ao = 0. Then, 
provided the eigenvalues of A all have negative real parts, we can show that the process 
{Xt} defined by 

Xt^<j f e^(*-"),5pdi?f (6) 

is a strictly stationary solution of (3) for t G (— oo, oo) with the corresponding CARFIMA 
process 

Yt=a f /3'e^(*-")5pdi?f. (7) 

The proof of (6) is similar to that of (4) and is hence omitted. 

For 1/2 < H < 1 and /,5 e L^{R;R) n L^{R;R), Gripcnbcrg and Norros [20] proved 
that 



/ pQC POO 

cov / /(«)di?f, / g{v)dB^ 

\'J — OO 'J — OO 



/ f{u)g{v)\u-vf 

-OO ^ — OO 



(8) 



For H < 1/2, Norros et al. [27] defines a class of stochastic integrals with the fractional 
Brownian integrator and deterministic integrands that are functions of bounded variation. 
However, they only consider integrals over the interval [0,cxd). In our case, integrals 
over an interval of the form (— oo,t] for finite t arc needed. For the definition of such 
stochastic integrals, we outline an approach similar to Norros et al. [27]. Let t be a 
fixed but arbitrary finite real number. Let F denote the integral operator mapping a 
boimded- variation function /(s), s <t, to a hmction r/(s), s <t, defined by the equation 

rf{s)^Hf{t)\t-s\'"-^sgn{t-s) 

+H f \u-s\^"~^sgn{s-u)df{u), 



where sgn(u) equals —1, or 1, depending on whether u is negative, zero or positive. 
Next, define an inner product between two functions / and g by the formula 

{f,9)r = / g{v)Tf{v)dv. 
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Define Lp(— cx),i] as the space of equivalence classes of measurable, bounded- variation 
functions / such that (/, /)r < oo. Consider the association of the simple function 
l[a,b] to — that preserves the inner product, which follows from equation (9) 
below, and hence the association can be extended to an isometry between the Gaus- 
sian space spanned by B^ ,u < t, and the function space ip(— oo,t] so that the integral 
jloo •^(") '-^^u '^^^^ defined as the image of / in the isometry. Consequently, provided 
that the integrals of the right-hand side of equation (9) exist, then for < H < 1/2 and 
all s, t£ R, 



cov(r /HdBf, / giv)dBf 

^Hf{s) [ \s~v\^"-^sg-o.{s~v)g{v)dv 



— OO 



H / / g{v)\u-v\^"-'-sgn{v^u)df{u)dv. (9) 



To demonstrate the above formula, it suffices to consider the case t^ s. The proof of 
equation (9) for simple functions is given in Section 3 and its validity for general functions 
then follows from the isometry alluded to above. 

From the above discussion, we can derive the stationarity condition and the autoco- 
variancc function of the CARFIMA process with Q < H <1. We state the stationarity 
condition of the CARFIMA process in Theorem 1. 

Theorem 1. Let < H < I. Equation (4) with a deterministic initial condition admits 
an asymptotically stationary solution if and only if all the eigenvalues of A have negative 
real parts. Moreover, under the preceding eigenvalue condition of A and assuming the 
solution is stationary, Yq and {B^ ,t > 0} are jointly Gaussian with the covariances 
given by 

/>oo 

cov(yo,Sf) = ifa/3' / c^"Sp{{u + t)^"-^ ^u^"~^}du, (10) 
Jo 



and the stationary mean of {Yt} equalling fiy ~ ^olqI ci 



For a random initial condition Yq , which may be correlated with the fractional Brown- 
ian innovation process, it can be verified that the sufficiency part of Theorem 1 continues 
to hold if Yq lias finite variance. Furthermore, the theorem implies that, under stationar- 
ity, yo a-nd the fractional Brownian innovations B^ ^ i > 0, arc generally correlated when 
B 7^ 1/2. This contrasts with the case B = 1/2, where the stationary distribution of Fq 
is independent of the standard Brownian motion. 

In Theorem 2, we use the covariance formulae of stochastic integrals given in ex- 
pressions (8) and (9) to calculate the autocovariance function of {Yt). In part (a), the 
autocovariancc function of the CARFIMA process is expressed in terms of three integrals. 
In part (b), the eigenvalues of the companion matrix A are distinct, so we have a closed 
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form of the autocovariance function. Then, (c) describes the asymptotic expression for 
the autocovariance function with H ^ 1/2. 

Theorem 2. LetQ<H <l. 

(a) Under stationarity, for h>0, the autocovariance function of {Yt} equals 

-fY{h) := cov {Yt+h,Yt) 

= H(3'Ae'^^(^J^ e~'^''u^"-Uu^V*P (11) 

-Hf3'Ac-'"'( r c''''u^"-'du]v*l3~Hl3'Ac^''( /°° c^V^-Mu ) 



where V* = e'^^'SpS'pB^'^ du. 

(b) Under stationarity, and when the eigenvalues Ai,...,Ap of the companion matrix 
A are distinct, 



ih) = ^r{2H + 1)^2 f/wl^f A., h), (12) 



where h>0, a{z) = — a^z^ ^ — • • • — ai, a^^\-) denotes its first derivative, /?(z) = 

1+ PlZ + I32Z'^ + ■ ■ ■ + PqZ'i , 

u{H, A, h) = 2{-Xf-^" cosh(A/i) + \^-^" c^^ P{2H, Xh) 

(13) 

and P{a,z)^ e^^u'^^^ du/T{a) , where the integration is along the radial line in the 
complex plane from to z. 

(c) For H ^ 1/2, as h ^ oo, we have the asymptotic expression 



jyih)r.a^H{2H-l)^h'"-^ (14) 

where "~ " means that the ratio of the left- and right-hand sides converges to 1. 
For H = 1/2, it is easy to verify that Theorem 2(a) and (b) can be simphfied to 

/3(A.)/3(-A.) 
^ ^a(i)(A.Oa(-A.)' ' 

which is consistent with the autocovariance function of the short-memory CARMA model 
proposed by BrockwcU [11]. Let r(a, z) = e~^u'^~^ du be an incomplete Gamma func- 
tion with complex arguments. Then, the fact that T{a -\-l,z) = aT{a, z) -\- z°'e~^ can be 
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used to show that, for 1/2 < H < 1, equatfon (13) is the same as expression (6.4) of 
Brockwell and Marquardt [12] up to a factor that involves the Hurst parameter only. 
The asymptotic expression (14) implies that for < < 1/2, 7y(/i) < when /i — *■ oo, 
which shows that the CARFIMA(p, ff, g) model with Q < H < 1/2 is antipersistent. In 
contrast, the model with 1/2 < H < 1 is of long-memory type. 

The autocovariance function established in Theorem 2 can be used to compute the 
spectral density function of {Yt,t G R} stated in Theorem 3 below. 

Theorem 3. For < H < I, the spectral density function of {Yj, t G R} is given by 

/y(«;)==^r(2i/+l)sin(7ti/)|«;|i-2^^^, (-00,00). (15) 

For < H < 1/2, equation (15) implies that /^7F(T)dr = 27T/y(0) = 0. The equa- 
tion also shows that the spectral density function of the CARFIMA(p, H, q) process is 
essentially a product of the spectral density of the ARMA process and the spectral den- 
sity of the fractional Gaussian noise. Thus, the CARFIMA(p, H, q) model is generated by 
applying an ARMA filter to the fractional Gaussian noise. Furthermore, the CARFIMA 
model is antipersistent if the fractional Gaussian noise is antipersistent, whereas it is 
long-memory if the noise is long-memory. Compared to fractional Gaussian noise, the 
CARFIMA model displays a much wider spectrum of autocovariance patterns, including 
non- monotone autocovariance functions. 

One major problem with continuous-time modeling is the identifiability of the 
continuous-time model, given discrete-time data. Let {Yih}i=i....^N be the observations 
sampled from a stationary C ARFIM A{p, H,q) process, where h is the step size. By the 
aliasing formula (Priestley [30]), the spectral density of {Fi/i}i=i,...,Ar equals 

f,iu;;6y) = lY.fy(^^^^)' -G[-7t,7t], (16) 
kez ^ ^ 

where /y(-) is as defined in equation (15). Using the frequency domain method, Tsai and 
Chan [38] showed that the CARFIMA (p, 7?, g) model with 1/2 < < 1 is identifiable, 
given regularly spaced discrete-time data. Specifically, they showed that, for (Oi^af) ^ 
(02,0-2), the set {uj\fh{w;Oi,af) ^ //i(a;; 02,^2)} positive Lebesgue measure if 1/2 < 
H <1. Because the spectral density function of {Yih}i=i,...^N is given by the same form 
as (16) for < H < 1/2 and 1/2 < 7? < 1, it can be shown by similar arguments that the 
CARFIMA model with < < 1/2 is also identifiable (see Section 3). The identifiability 
problem with H = 1/2 is more difficult; sec [38] for further discussion. In summary, we 
have the following theorem on identifiability. 

Theorem 4. Let Y = {Yt.}fLi be sampled from a stationary (Gaussian) CARFIMA (p, 
H, q) process given by equation (1), where < H < 1, H ^ 1/2, a{-) and /?(•) have no 
common zeros, all roots of a{z) ~ and the roots of P{z) ~ have negative real parts. 
If the step size ti ~ ih with h > 0, then for {61, af) ^ (02,o'|), the set {uj\ffi{uj;0i,(jf) 7^ 
/ft,(a;; 6*2, erf )} has positive Lebesgue measure. 



186 



H. Tsai 



We note that the roots of the polynomial a(-) are the same as the eigenvalues of the 
matrix A and the condition on the roots of a{z) = is necessary for the stationarity of 
the process, whereas the condition on (3{z) = is akin to the invertibility condition for 
discrete-time processes. 

3. Proofs 



Proof of equation (9). Consider simple functions of the form /(m) = X^^g Cil(s. s-^^j (u), 

where Sm = s, and g{v) = X]j=o^ '^jl(tj,t3+i](^)j where t„ = t. Then, the left-hand side of 
(9) becomes 



Covf r f{u)dB^,f g{v)dB^ 

\J ~oo J —oo 



(rn—l n — 1 \ 

^ c.«^^ - < ), d,{B^^^^ - B«) (17) 
i=0 j"=0 / 



ni— 1 n— 1 

2 



i=0 j=0 



The first term of the right-hand side of (9) is 

J — OO 



"Z^l rtj + i 

= HCrn-l22'^i \Sni-v\'^" ^Sgn{Sm-v)dv (18) 

-i—n '^ii 



3=0 
n-1 



3=0 

If we let c_i =0, then the second term of the right-hand side of (9) is 
H f [ g[v)\u-v\'^"~^sgn[v-u)df{u)dv 

J —oo J —OQ 

/t ™-l 
giv) E ~ v\^"~'^ sgn(t) - Si){ci - Q_i) di; 

m-l n-1 (-tj + i 

= i/E(C'-C'-i)E^^_/ \s^~v\''"-Hgn{v~s,)dv (19) 
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m—l n — 1 

= i ^ (C, - C,_i) dJ{\s^ - tj + l\'" - |S. - t,\'"} 
1=0 j=0 



m—l n — 1 

2 



i=0 j=0 



i=0 j=0 

Therefore, by (18) and (19), the right-hand side of (9) becomes 



m—l n— 1 

2 



=0 j=0 



which is the same as equation (17). This proves the vahdity of equation (9) for simple 
functions. □ 



Proof of Theorem 1. The proof for = 1/2 is triviaL For the proof where 1/2 < H < 1, 
see Tsai and Chan [37]. We now consider the case where < H < 1/2. The proof of the 
first part of the theorem is similar to that of Theorem 1(a) in [37] and is hence omitted. For 
the proof of equation (10), first note that equations (4) and (6) are essentially equivalent 
if ao ~ 0. Since, by (6), we can write Xq ~ a /"^ e~^"(5p dS^ and Bt = J^dB^ , by 
equation (9), we have 

cov(yo,sf) 

= cov (^a £ dSf , £ dSf ) 

= Ha(3'Sp f \-v\^"-'^sgn{-v)dv (20) 
Jo 

-Ha f f \u~-v\^"-^sgn{v-u)(3'Ae-'^''5pdudv 

Jo J-oo 



10 J-oo 

noc 



-I3'5j,t'" - - / {{t + uY'' -u'"}P'Ae''^5pdu. 







Now, based on the integration by parts technique, the equality in (10) follows from 
equation (20). The stationary mean of {Yt} follows from equation (5) and the subsequent 
discussion. □ 
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Proof of Theorem 2. The proofs of 2(a) and (b) with H = 1/2 are trivial. First, we 
prove part (a) with < < 1/2. By equations (7), (9) and routine calculus, we have 

7y(/i) = a^H f /3'5p(5;e'^'(*-")/3(i + h- vf"-^ dv 

J — OO 

/t pOO 
/ (3'Ae^^*+''+^-^'^SpP'e''^'-^'^dpw''"-'dwdv (21) 
^oo Jo 

/t pt-\-h—v 
/ p' Ae^^'+'^-^-^^SpP' c^^'-^^dpw^"-' dwdv. 

-OO Jq 

Now, equation (11) follows from equation (21) and equation (6.20) of Karatzas and 
Shreve [22], namely, AV* + V* A' = -cr^SpS'^. For l/2< H < 1, equation (11) follows 
from Theorem 1(c) of Tsai and Chan [37] and the integration by parts technique. 

Below, we prove part 2(b) with < H < 1/2. The proof with 1/2 < iJ < 1 is similar 
and is hence omitted. When all the eigenvalues of A have negative real parts and are all 
distinct, Brockwell and Marquardt [12], equation (2.15), show that, for h>0, 

Differentiating the above equation with respect to h on both sides, we have 

f^'^-'% = t§^/'''- (23) 
By expressions (21), (22), (23) and routine calculus, we obtain 



7y(/i) 



(24) 



--r(2i/+l) 



.-i—l ^ ' 



^(1 



/3(A,) 



^^a(i)(A,)(A, + A,) 



Y-^A,) 
2-^ ^(1) 



(25) 



x{e-^^''(~A,) 
p 



X 



51 ^(1 



/3(A, 



-a(i)(AO(A, + A,) 
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By the partial fraction (FeUer [18], page 276) 

/3(A,-) 

«(«) ,tt«^'HA,)(.-A,)' ^ ^ 

we have 

/?(-A0 _ /?(Aj-) 

«(-A,) ^tt"^'HA,)(A. + A,)- ^ ^ 

Equation (12) now follows from equations (25) and (27). This proves (b). 

To prove (c) with < H < 1/2, we note that P{a, z) = 1 — r(a, z)/r(a). Furthermore, 
by equation (6.5.32) in [1], T{a,z) ~ z°~^e~^{l + (a — l)/z} when z — > cx). Thus, as 
h ^ oo, we have 

Hence, by equations (12) and (28), a.s h ^ oo, wc have 

jyih) ^ -2a-H{2H l)h^"-- (nix f ' ^^9) 

^ a^^i[Xi)a{-\i)\i 

Now, = l/{Ai(Ai + Aj)} + l/{Aj(Ai + A^)} and expression (27) implies that 



p 

E 



-«W(AOa(-A,)A, 



EE 



/3(A,)/3(A, 



-^«(i)(AOa(i)(A,)A,(A. + A,) 



1=1 ] 
p p 



1 - - /3(A.)/3(A,) 
2fe,tt«^'HAi)aW(A,)A.(A. + A,) 

ly-y /^(A.)/3(A,) 



(30) 



p p 



/3(A.)/3(A,) 



1 ^ ^ 

"2lt§«W(A.)aW(A,)A.A, 

1 /3^(0) 
2a2(0)' 
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where the last equality follows from the identity (26). The proof of (c) with < i? < 1/2 
now follows from expressions (29) and (31). Note that the proof of (c) with 1/2 < H <1 
is similar to that of0<i/<l/2 and is therefore omitted. □ 

Proof of Theorem 3. For the proof with 1/2 < H < 1, readers can refer to Tsai and 
Chan [38]. For case where H = 1/2 is trivial, the proof with G<i7<l/2 is similar to 
that of Theorem 1 in [38] and is hence omitted. □ 

Proof of Theorem 4. For the proof with 1/2 < H < 1, see Tsai and Chan [38]. The 
following proof with 0<-ff<l/2is similar to that in [38]. 

If we write fc{uj) = fyii^) ^ g{H)L{uj)\uj\^-^" , where g(i/) = r(2i/ + 1) sin(7TiJ)/(27t), 
L{lo) = a2|/3(ia;)|V|a(itj)|2, then fi^\uj) = g{H)\uj\-^" {\lo\L'^^\uj) + (1 - 2H)L{uj)} and 



log(/(i)(^)) = \og{f['H^)) + log 1 + ^ 



(i)( 

/ 

= \ogg{H) - 2ff log|w| + \og{\uj\L^^\uj) + (1 - 2H)L{uj)} (31) 
+ log{l + R{lu)), 

where i?(a;) = i?2(w)/i?i(a;), Ri{lo) = f^^\uj) and i?2(w) = 'Ek^o fc^\'^ + 2kn). It follows 
from lim;j^o{log(/^^-'(w))/log|ti'|} = —2H that the Hurst parameter H is identifiable. 
Following Tsai and Chan [38], it suffices to show that, given f{uj),uj G [— tt, tt], we can 
determine L and all of its higher derivatives at w = 0. Below, we show that L{0) and 
L(i)(0) are identifiable. The identifiability of i('=)(0) for fc > 2 can be proven similarly by 
the arguments in [38]. The identifiability of L(G) simply follows from the fact that 

lim {log(/(i) (lu)) - loggiH) + 2H\og |co| - log(l - 2H)} 



: lim 



-logi(O). 

Next, we prove the identifiability of i?2(0), which is needed to prove the identifiability of 
L^^\Q). By equation (31), we have 

^{log(/(i)(c.)) + 2if logical} 



\u\L^'^\lu) + [2 - 2H)L(^\uj) R^^\uj) 



where 



|a;|L(i)(w) + (l-2i/)i(w) 1 + R{ujY 



+ R{u;) (l + i?(^))i?i(c.) {r,{cu)r['\lu) R^^\uj) 



(32) 



(33) 
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lim 



\1~2H 



duj 



{log(/(i)(^)) + 2i/log|^|} 



lim 



uj^o[ 1 + R{lu) 



- lim 



lim 



Rliuj) 

R2iLj){\Lj\^L^^^uj) + 2(1 - 2H)\uj\L'-^\u;) - 2i7(l - 2H)L{uo)} 



g[H){\u;\Li^){Lo) + [l~2H)L{u,)y 



2HR2{0) 



g{H){l-2H)L{Qy 
hence -R2(0) is identifiable. Now, equation (33) implies that 
i?(i)(tj) 2HR2{0)\uj\^"-^ \ 



ci™ [ 1 + R{uj) g{H){\ - 2H)L{Q) j 

R2{uj)Rf\u) 2HR2mu:\^"-^ \ 
Rliuo) g{H){l-2H)L{0)\ 



■ lim 



: lim 



2H 



R2{Lu)\uj\'^"~^{\io\'^L^'^\uj) + 2(1 - 2H)\u\L''^\uj) - 2H{1 - 2H)L{uj)} 
g{H){\u;\L^^){u;) + il-2H)L{cu)}^ 

2i/i?2(0)|c^|2»-i T 



(34) 



giH)il-2H)L{0) 



g{H) u>^a 




\2H-1 



^2(0) ] 

{\uj\LW{uj) + (1 - 2i7)L(a))}2 (1 - 2H)L{Q) J 



(l-2F)L(w)i?2(w) 



because i?2(^^) = i?2(0) + wO(l) and i(w) = L(0) + a;0(l) for w tending to 0. Therefore, 
by (32), (33) and (34), 



lim 



A{iog(/(^)(.)) + 2i^iog|.|} - 2/^^^(0)1-1^;; ' 

j.^^^ |a'|L(2)(u;) + (2-2g)L(^)(u;) 
"t^™ |cj|L(i)(w) + (1 -2iJ)L(tj) 

(2-2ij')£(^)(0) 
{1-2H)L{0) ' 
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which proves the identifiability of L(i)(0). □ 

4. Conclusion 

We have proposed a unified continuous-time framework that is useful for studying time 
series with short memory, long memory and antipersistence. The identifiability of the 
CARFIMA process with discrete-time data established in Theorem 4 is a fundamental 
feature that makes the model practical for data analysis. Therefore, in future research, 
it would be both interesting and useful to study the statistical inference of antipersistent 
CARFIMA models for data analysis. 
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